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Statement of Work 


Object ive 

Tl'i*? objective of this investigation is to evaluate the applicability of 
el(>ctromagnctic deep-sounding experiments using natural sources in the 
magnetosphere by incorporating Magaat data w’th other geophysical data. 

Approach 

The Investigator shall pursue the above objective through an analysis of 
Magsat satellite data, ground-based magnetic observations, appropriate 
reference field models, and other satellite data. 

The objective will be pursued by seeking the optimal combination of 
observations v;hich lead first to a global, and then to a regional, 
character izatlon of the conductivity of the Karth's upper u'antle. 

Tasks 


The following tasks shall be performed by the investigator in fulfillment 
of the above objective: 

a. Use data from Magsat satellite to constrain a long-period 
global "response function" for the average Earth at low latitudes over 
a period ranging from C hours to 27 days. 

b. Synchronize the Magsat data with low-latitude ground-based 
observatory data to determine the vertical gradient of the respective 
magnetic field components. Use the vertical gradient of the appropriate 
components to Independently ascertain the separation of external and 
internal field contributions. 

c. Segregate tl:e Nagoat electromagnetic "response functions" 
according to the tectonic regime at the Earth's surface and evaluate 
systematic differences between regions having lateral scale sizes on 
the order of 1000 km or greater. 

d. Theoretically evaluate problems of resolution and interpretation 
involving electromagnetic induction by temporally and spatially-varying 
magnetospheric sources in a rotating Inhomogeneous Earth as observed at 
arbitrary points in space. Use these tlieoretical studies to constrain 
the interpretation of Magsat data as well as to propose further applica- 
tions of satellite-based electromagnetic deep-sounding experiments. 

e. Integrate the regional response functions with other geophysical 
data in order to constrain the joint interpretation of comprehensive 
physical models. 

f. Prepare and submit to NASA periodic progress reports and a 
detailed final report documenting the results of this investigation. 
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Our preliminary goal in the MAGSAT project is to coordinate ground-based 
and satellite data sets. Toward this end we have developed, and are in the final 
stages of testing, a spherical harmonic analysis program which takes magnetic 
data in universal time from a set of arbitrarily spaced observatories and cal- 
culates ii value for the instantaneous magnetic field at any point on the globe. 
The calculation is done as a 1<. ast mean-squared value fit to a set of spherical 
harmonics up to any desired order n. 

The program is also designed to accept as a set of input parameters the 
orbit position of a satellite and to coordinate it with ground-based magnetic 
data for a given time. The output is a predicted time series for the magnetic 
field on the earth's surface at the (r,0) position directly under the hypothet- 
ical orbiting satel Ite foi the duration of the tli. 3 period of the input data 
set. 

Using this program to "track" the vsurface magnetic field beneath the sat- 
ellite will allow one to compute narrow-band averaged crosspowers between the 
spatially coordinated satellite and the ground-based data sets. These cross- 
powers can then be used to calculate field transfer coefficients with minimum 
noise distortion. As an example, we shall discuss the application of this 
technique to calculating the vector response function, W. 

The following variables represent a narrow-band filtered frequency 
spectrum data set of the respective satellite magnetic component indicated 
(we assume the NEV coordinate system); 

HgAT* north component of satellite data 
ZgAT* vertical component of satellite data 

Next, if we assume a external field, we can define: 
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Z “ ZeA.i./c 08 0 

C hAi 


U) 


to be the ' latitude-compensated" compoiientH of the satellite data. 

After Banks (1969), we define a vector response function W, where: 


Zc “ W»c <2) 

The usefulness of W centers on Its Independence of source field details. 
This result cun be obtained from multi-layered spherical earth models, as long 
us the source field assumption is maintained. 

Noting the form (2), we introduce the term W', a least-squares estimator 
of W; if we let V stand for the statistical variance of any estimation of W, 
say W , from W, then the following equation stands as our definition of W' : 

dV/dW - - W'H^|^>), 

“ 0, (3) 


where the angular brackets denote a frequency band average of the form: 


<(Z -W'H )(Z -W'H )*> 
c c c c 


J (lo) -W ' (w) ) (Z^ (w) -W ' 11^ (to) ) *d (0 

(o^lAto/Z 


where 


(OqI center frequency of the respective filter band, 

Auk the band width of tlie filter, 

Z^(u)), H^(ui): elements of the respective spectra, Z^ and 

In general we shall term "<AB*> the "inner product" or "crossproduct" 


of A and B. 
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As a result of (3) we find that: 


<Z H*> 

<H H*> 
c c 


(4) 


In a noise-free situation, W* would represent 

up to the resolution of our instruments. \p it is, 

contain noise, such that H and Z may be expressed 

c c 


an exact calculation of W, 
however, our data will 
as 


H - + N„ 

c c H 

Z - Z**^® + N,. 
C c Z 


(5) 


where and are tlie noise spectra associated with the M and Z components, 
respectively. 

We may evaluate how (4) will be affected by the noise. First, we see that 
the numerator, that is, the inner product of Z and H, will be essentially unaffec- 
ted. We compute: 

<Z H*> =» <(Z**^® + N„)((H‘^^®)* + N;'s)> 
c c ^ c H 

- <7/^®!!!'^^^®)* + N„N* + , ‘''^®N* + N„(ll^^®)*> 
c c Z H c H Z' c 

<7®^®(h”^®)a> (6) 


since approaches zero because of the random and uncorrelated nature of 

noise sources. For the same reason, <z‘^^®N{;> and <N„ (H"'*^®)*> approach zero as 
well, so long as the selectivity of the band-limiting filter applied to the 
data is wide enough to allow the noise to average out. 

Therefore, expression (6) supports our assertion that the numerator of 
(4) is unaffected by noise. The denominator, on the other hand, is decidedly 


affected: 
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<iyi*> - <(nf K + -f 


<hif^!* + |N„|^ 


( 7 ) 


where the cross terms drop out accordlnp, to the same reasoning employed In cal- 
culating (6). 

The total effect of noise on o..r estimated parameter W' In (4), then, is 
a downward bias in magnitude (the denominator is real) relative to its noise- 
free value W. In other words, if we say 


^^siy jSig)^^ 
c c 

c c 


(8a) 


then 


W 




<{H 


sl8 


c _ _ _ 


.'i w. 


(8b) 


the equality holding in the limit as N^j approaclies zero. 

We note that only the noise in the H data affects the result. In effect, 

c 

in formulating (2), we cast in the role of the "Independent parameter", from 
which could be obtained by way of tl>e transfer function W. This was accept- 
able e.\cept that the noise in H, unaccountable fo in our analytical transfer 
function, biased our final estimate. 

Let us try a simple recasting of roles whereby "becomes" independent 
and noise-free in the foxmulation. We define Y such t!iat 
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where In a nolBc-frec Bltuation, 

1/Y - Z' - E. 

Performing a least-squares fit for Y as we did above for W* , we obtain: 
<H Z*> 


V - 


c c 


c c 


(9a) 


Then, performing a signal-noise separation as we did boding to (6) and 
(7) , we find: 


Y» 


+ 1nJ2^ 


-C • • 'Z' 

This result is analogous to (8b); the magnitude of Y' will be biased 
downward by the noise in Z , N,. , as W was biased downward by N„ . If we let 


W" - 1/Y' 

W" will then be biased upv/ard by C!ic N„. Hence we have: 

4i 

w'<w<w" 

as a constraint on our hypothetical noise-free estimate, W. 

A possible improvement would be to take the geometric mean of W and 
W", from (8b) and (9b) respectively, resulting in: 


( 10 ) 


( 11 ) 


2 (W'W")^'^^ - ( ) • ( 






+ N. 




) 


(12a) 


If we express; 


c ' c ‘ 

„Sig , ij^Sig, i({)M 
c ‘ c ' 

then using this notation and the reasoning leading to (6), we have: 


<|z®^S|=' + In 


-) 




N„l> 

Examination shows that is an excellent for the the magnitude of W if 

InJ' |n„M 


^sig I 2 




( 13 ) 
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Ill addlC:Lon, phase information is worried noise-free in w , and is 
shown to be« equivalent to, the difference in the phase angles of the 
"hypothetical" noise-free vali.es, and Caution must be 

maintained, however, in using W^y for dato gathered at high colatitudcs, 
where is in general much smaller than In this ease the 

approximate equality in (13) will break down. The necessary correction 
would involve reformulating (l-^a) in the form of a weighted 

geometric mean, where the weight Is a function of the appropriate coherencies. 


The above discussion considered a metiuni for extracting W from a data set 
from one site (a satellite in our Instance). As we mentioned at the beginning 
of this discussion, the addition of a data set from another site, coordinated 
with the primary data set, will afford us a far more reliable method for extract- 
ing noise-free parameters and, in particular, fi>r formulating W. Let us define 

* H'''";/sin 0 
ref ref 


Analogous to our definitions in (1), represents the latitude-corrected 

narrow-band frequency spectrum for the horizontal, northern component of a 
ground-based data set, spatially and secularly coordinated with our satellite 
components and 

Let us define a new field parameter as follows: 


X 


<Z ri* 
c ref 

<H H* .> 
c ref 


(15) 


where we have teken (4) and replaced ir< with 


«?.£ ■ <CP* ^ Mr- 


We make the separation 


and assume that there Is no statistical correlation between ground-based noise 
and satellite-based noise. It follows then from our arguments leading to (6) 


that 



7 


, c ref 

c r^f 


(16) 


Hence, (15) mlnlmizeb the cfiectH of noise in an csClmato of X. However, wc 

must address the question of the physical interpretation of X; can it be equated 

with W? More specifically, what assumptions must we make about the relationship 

between and so that "X “ W" nuiy be asserted? Another way of posing the 

question would be; in a noise-free situation, what constraints must be placed 

on H'*' c so that its substltutlc-n for H* into (4) would not disturb the identity? 
ref c 

Mathematically, the answer is simple: the frequency spectrum represented 
by must be a scalar multiple of the frequency spectrum represented by 

for the substitution to be valid. In physical terms, this leads us to two basic 
constraints; one relates to the properties of the space between the ground and 
th«? satellite and the other relates to the conditioning of the data set itself. 

For a spherically symmetric layered earth, assuming a source field, 
it can be shown that at or above the uppermost boundary, 

- 11 + |(r^/r)^j™ (17) 


where r^^ is the radius of the uppermost boundary, r is the radius of the observer 
(it is presumed that r > r^^), is the magnitude of the external driving 
field, and R is a complex response function: 

(iuprj^ - 27 /) 

~rr 


R 


(iwprj^ + Z') 


(18) 


Z' is the surface impedance, calculated by iteration from the bottommost surface 
of the model. 

The space between ground and satellite must have no sharp media boundaries 
and the conductivity throughout must be near-zero. Or else, reflection and/or 
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attenuation of tue field betweon ground and •^atullltc will occur In an uneven 
nuinner across the band oT the frequency specira of the magnetic components. 

The data set conditioning must have a selectivity sufficiently narrow to 
insure that the response of the earth ns a whole is essentially constant over 
the filter bandwidth. This last constraint i-v iresents a tradesff with the 
assumptions made to reach (6), (7), (12) and (16), our various noise-minimized 
"W-estimators" . We therefore must bo aware of r'.;,* special care necessitated in 
the choice of a proper selectivity, especially with long-period data. 

We have presented these constraints under which "X ■ W" is a good approx- 
imation in a very qualitative manner. However, they are certainly basic problem* 
that must be taken into account when applying any quantitative model to the data. 

An attractive feature of the formulation proposed in (16) Is the absence 
of a vertical ground-based component. This is an extremely useful property of 
this analysis as it is well known that surficlal lateral inhomogeneltles will 
introduce a distortion in the surface vertical component of far greater 
magnitude than the distortion introduced in the liorlzontal components. 

In conclusion, we make two observations. First, as W is Independent of 
source field details, its value as calculated from (16) should remain relatively 
constant as different segments of the 7-odd months of MAGSAT are processed via 
this formula. The size of the deviation of a set of "W's" thus calculated will 
provide us with a quasl-quantitative measure of the validity of the source 
field assumption which underlies the vector response function's source-independent 
character. This assumption can also be checked with ground-based data. 

Finally, we note that although long-period data is of great use in our 
analysis as it suffers least from attenuation in the atmosphere or distortion 
by surficial inhomogeneities, it is also most difficult to extricate from the 
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dii'Ca baac. Saccllite data dost! to tho auroral zone nuist bo discarded because 
of the large offo.c of the field-aligned currents have on the data. Therefore, 
we can have at most about 100 degrees of orbit time (about 23 minutes of data) 
for a continuous data sot. We ire at this time still looking tor alternative 
methods of chaining data sets together whicii minimize possible effects of 
spatial aliasing. 


